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ABSTRACT 

We study the baryon mass fraction in a set of hydrodynamical simulations of galaxy clus- 
ters performed using the Tree+SPH code GADGET-2. We investigate the dependence of the 
baryon fraction upon the radiative cooling, star formation, feedback through galactic winds, 
conduction and redshift. Both the cold stellar component and the hot X-ray emitting gas have 
narrow distributions that, at large cluster-centric distances r > -R500, are nearly independent 
of the physics included in the simulations. Only the non-radiative runs reproduce the gas frac- 
tion inferred from observations of the inner regions (r w -R2500) of massive clusters. When 
cooling is turned on, the excess star formation is mitigated by the action of galactic winds, but 
yet not by the amount required by observational data. The baryon fraction within a fixed over- 
density increases slightly with redshift, independent of the physical processes involved in the 
accumulation of baryons in the cluster potential well. In runs with cooling and feedback, the 
increase in baryons is associated with a larger stellar mass fraction that arises at high redshift 
as a consequence of more efficient gas cooling. For the same reason, the gas fraction appears 
less concentrated at higher redshift. We discuss the possible cosmological implications of our 
results and find that two assumptions generally adopted, (1) mean value of Yj-, = n J^ n not 
evolving with redshift, and (2) a fixed ratio between / star and / gas independent of radius and 
redshift, might not be valid. In the estimate of the cosmic matter density parameter, this im- 
plies some systematic effects of the order of Af2 m /f2 m < +0.15 for non-radiative runs and 
Af2 m /f2 m ss +0.05 and < —0.05 for radiative simulations. 

Key words: cosmology: miscellaneous - methods: numerical - galaxies: cluster: general - 
X-ray: galaxies. 



1 INTRODUCTION 

Baryons in galaxy clusters, mainly in the form of stars in galaxies 
and hot X-ray emitting plasma, trace the potential of the collapsed 
structure in which they fall into. At the same time, their spatial dis- 
tribution and thermodynamical properties are affected by the phys- 
ical processes acting on them. Measurement of direct cluster ob- 
servables, such as stellar optical light and hot gas X-ray emission, 
is the only viable approach to investigate the physics that drives the 
evolution of such structures. The observed relations among proxies 
of the internal energy and entropy levels show that, at least in mas- 
sive systems, the dominant baryonic component of the intra-cluster 
medium (ICM) has a luminosity, temperature and mass that well 
follow those predicted under the assumptions of a plasma emitting 
by bremsstrahlung and sitting in hydrostatic equilibrium with the 



underlying dark matter potential (e.g. Kaiser 1991, Evrard & Henry 
1991). 

Moreover, the fact that galaxy clusters are relatively well iso- 
lated structures, that form in correspondence of the highest peaks 
of the primordial gravitational fluctuations, suggests that their rela- 
tive baryon budget and the mass function are highly sensitive tests 
of the geometry and matter content of the Universe. In particular, 
using X-ray observations of the baryonic content to infer the gas 
and total mass of relaxed clusters in hydrostatic equilibrium allows 
one to place a lower limit to the cluster baryon mass fraction, which 
is expected to match the cosmic value fib/Hm (White et al. 1993, 
Evrard 1997, Mohr, Mathiesen & Evrard 1999, Ettori & Fabian 
1999, Roussel, Sadat & Blanchard 2000, Allen et al. 2002, Ettori et 
al. 2003, Ettori 2003, Allen et al. 2004). 

The purpose of the present work is to use an extended set 
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of hydrodynamical simulations of galaxy clusters, treating a va- 
riety of physical processes, to study how the spatial distribution of 
the baryons, as contributed both by the stellar component and by 
the hot X-ray emitting gas, are affected by the physical conditions 
within clusters. In this respect, our analysis extends previous anal- 
yses of the baryon fraction in cluster simulations which included 
only non-radiative physics (Evrard 1990, Metzler & Evrard 1994, 
Navarro et al. 1995, Lubin et al. 1996, Eke et al. 1998, Frenk et 
al. 1999, Mohr et al. 1999, Bialek et al. 2001), and the processes of 
cooling and star formation (Muanwong et al. 2002, Kay et al. 2004, 
Ettori et al. 2004, Kravtsov et al. 2005). 

The paper is organized as follows: in Section 2, we describe 
our dataset of simulated galaxy clusters; in Section 3, we present 
the results obtained from this set of simulations on the gas and 
baryon fraction, also comparing these results to the observational 
constraints. In Section 4 we compare our results to previous analy- 
ses based on different simulations. Finally, we summarize and dis- 
cuss our findings in Section 5. 



2 PROPERTIES OF THE SIMULATED CLUSTERS 

We consider two sets of clusters, which have been selected from 
different parent cosmological boxes. The first set is extracted from 
the large-scale cosmological simulation presented in Borgani et al. 
(2004). The second one is a re-simulation of 9 galaxy clusters, ex- 
tracted from a pre-existing lower-resolution DM-only simulation. 
Our simulations were carried out with GADGET-2 (Springel 2005), 
a new version of the parallel Tree+SPH simulation code GADGET 
(Springel et al. 2001). It includes an entropy-conserving formula- 
tion of SPH (Springel & Hernquist 2002), radiative cooling, heat- 
ing by a UV background, and a treatment of star formation and 
feedback from galactic winds powered by supernova explosions 
(Springel & Hernquist 2003). 

2.1 Clusters extracted from a cosmological Box 

The first set of simulated clusters has been extracted from the large- 
scale cosmological simulation of a "concordance" ACDM model 
(Q m = 0.3, fi ( ,A = 0.7, ft b = 0.019ft -2 , ft = 0.7, cr 8 = 0.8; 
Borgani et al. 2004). Here we give only a short summary of its 
characteristics, and refer to that paper for more details. The run 
follows the evolution of 480 dark matter particles and an equal 
number of gas particles in a periodic cube of size 192ft" 1 Mpc. 
The mass of the gas particles is m gas = 6.89 x 10 8 ft _1 Af Q , while 
the the Plummer-equivalent force softening is set to 7.5ft -1 kpc 
from z = to z = 2, while kept fixed in comoving units at higher 
redshift. At z — we extract a set of 439 mass-selected clusters 
with virial masses Mvir > 5 x 10 13 ft _1 M o . 

The efficiency of conversion of the energy provided by SN 
explosions into a kinetic feedback (i.e. winds) is set to 50 per cent, 
which gives a wind speed of ps340 km s _1 . 

2.2 Resimulated clusters 

This set includes simulations of 4 high mass systems with 
M 2 oo =(l.l-1.8)xl0 15 ft _1 M©. The cluster regions were ex- 
tracted from a dark-matter only simulation with a box-size of 
479 ft -1 Mpc of the same cosmological model as the first set, but 
with a higher normalization of the power spectrum, ug = 0.9 (see 
Yoshida, Sheth & Diaferio 2001). Using the "Zoomed Initial Con- 
ditions" (ZIC) technique (Tormen 1997), they were re-simulated 



with higher mass and force resolution by populating their La- 
grangian volumes in the initial domain with more particles, while 
appropriately adding additional high-frequency modes. The initial 
particle distributions (before displacement) are of glass type (White 
1996). The mass resolution of the gas particles in these simulations 
m gllB — 1.7 x 10 8 ft -1 Mq. Thus, the clusters were resolved with 
between 2 x 10 6 and 4 x 10 6 particles, depending on their final 
mass. For all simulations, the gravitational softening length was 
kept fixed at e — 30.0 ft -1 kpc comoving (Plummer-equivalent), 
and was switched to a physical softening length of e = 5.0 ft -1 kpc 
at z = 5. 

For simulations including star formation and feedback the SN 
efficiency in powering galactic winds is set to 50 per cent, as in the 
cosmological box, which turns into a wind speed of ~340 km s~ . 
For the sake of comparison, some runs have also been performed 
also by switched off the winds completely or increasing the SN 
efficiency to unity. 

For some of the cluster simulations we also included the effect 
of heat conduction. Its implementation in SPH, which is both stable 
and manifestly conserves thermal energy even when individual and 
adaptive time-steps, has been described by Jubelgas, Springel & 
Dolag (2004). This implementation assumes an isotropic effective 
conductivity parameterized as a fixed fraction of the Spitzer rate, 
that we assume to be 1/3. It also accounts for saturation, which can 
become relevant in low-density gas. For more details on the proper- 
ties of simulated galaxy clusters including thermal conduction see 
Dolag et al. (2004). 

Furthermore, some simulations were carried out using a mod- 
ified artificial viscosity scheme suggested by Morris & Monaghan 
(1997), where every particle evolves its own viscosity parameter. 
Whereas in this scheme shocks are as well captured as in the stan- 
dard one, regions with no shocks do not suffer for a residual non- 
vanishing artificial viscosity. Therefore turbulence driven by fluid 
instabilities can be much better resolved and as a result of this, 
galaxy clusters simulated with this new scheme can build up suffi- 
cient level of turbulence-powered instabilities along the surfaces of 
the large scale velocity structure present in cosmological structure 
formation (Dolag et al. 2005). 

In summary, this set of cluster was simulated 4 times, with 
different kind of physics processes included: 

• Gravitational heating only (code=G). 

• Gravitational heating only with low viscosity scheme 
(code=GV): Like G, but using the alternative implementation of 
artificial viscosity. In this simulation, galaxy clusters are found to 
have up to 30 per cent of their thermal energy in the turbulent mo- 
tion of the ICM, leading to a sizeable contribution of non-thermal 
pressure support in the center of galaxy clusters; 

• Cooling + Star Formation + Feedback with weak winds 
(code=FwW): the wind speed is fixed at ^340 km s _1 . 

• Cooling + Star Formation + Feedback with weak winds and 
Conduction (code=FwWC): Conduction efficiency set to be 1/3 of 
the Spitzer rate. 

In order to have under control the effect of changing the feed- 
back efficiency, one high-mass cluster was further simulated with 
the following setups: 

• Cooling + Star Formation + Feedback with no winds 
(code=F): Like FwW, but with winds switched off. 

• Cooling + Star Formation + Feedback with strong winds 
(code=FsW): Like FwW, but with wind speed increased to «480 
km s _1 , corresponding to a SN efficiency of unity. 




Figure 1. The baryon and gas mass fractions at -R200 and z = as func- 
tions of the virial mass. The solid dots with error-bars refer to the mean and 
dispersion measured in the corresponding bins in virial mass. The number 
of objects in each bin is indicated. The most massive re-simulated systems 
are also shown. The shaded regions indicate the lcr range of Y" gas and Y\> 
for the set of simulated clusters extracted from the cosmological box. 

The center of each cluster is defined as the position of the par- 
ticle having the minimum value of the gravitational potential. Start- 
ing from this position, we run a spherical overdensity algorithm to 
find the radius Ra c encompassing a given overdensity A c , with 
respect to the critical one at the redshift under exam, and the mass 
Ma c enclosed within this radius. In the present work, we consider 
values of the overdensity A c equal to 2500, 500 and 200. The cor- 
responding radii relate to the virial radius, which defines a sphere 
with virial overdensity (of « 101 at z = and ~ 157 at z = 1 for 
our cosmological model and with respect to the critical value), as 
(^2500,^500,^200) ~ (0.2, 0.5, 0.7) x iivir- For each cluster, the 
hot gas mass fraction and the stellar mass fraction within a given ra- 
dius r are then calculated as / gas (< r) = M gas (< r)/M t ot(< r) 
and / st ar(< r) = M star (< r)/M to t(< r), respectively. 
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3 RESULTS FOR THE GAS AND STELLAR MASS 
FRACTIONS 

For the sake of clarity, we define the quantities Y gse , Y s t ar and Yt, as 
the ratios between / gas , / star and / b = / gas +/ 8 tai, and the cosmic 
value adopted in the present simulations, fib/fi m = 0.13. Their 
mean values (and standard deviations of the observed distributions, 
where available) at redshift z — 0, 0.3, 0.7 and 1 are quoted in 
TableQfor the different physical conditions considered. 

We compare in Fig.Qthe gas and baryon fraction as a function 
of the virial mass, for both the simulated clusters extracted from the 
cosmological box and for the subset of the re-simulated ones with 
M200 > 10 15 hT 1 Mq. This plot demonstrates that, when com- 
puted within the whole cluster virial region, there is no relevant 
dependence upon mass of such fractions. We remind the reader 



Figure 2. The radial distribution of the gas and baryon mass fractions (in 
unit of the cosmic baryonic value) at z = and 1 for one of the four mas- 
sive re-simulated systems. The vertical dashed lines indicate the location of 
^2500. ^?500 and -R200 f° r the Gravitational heating only case. 



that the two sets of clusters have been simulated for the same cos- 
mological model, the only difference being in the value of 0%, as- 
sumed to be 0.8 and 0.9 for the first and for the second set, re- 
spectively. In Fig.Q the re-simulated clusters indicated with a di- 
amond correspond to the same simulation physics as for the clus- 
ters extracted from the cosmological box. They show a mean Yb 
consistent with that measured for the massive systems of the cos- 
mological box, but with a slightly lower value of Y gas . This is due 
to the increase in the cooling efficiency with the amplitude of the 
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Figure 3. The gas (diamonds) and total baryon (dots) mass fractions, in unit of the cosmic baryonic value at R2500 (left panel), R500 (central panel) and -R200 
(right panel). For each physical case considered, we plot the mean and standard deviation values measured at z = 0, 0.7 and 1. The shaded region show the 
error- weighted mean and standard deviation of (1) / ga s (-R2500 ) estimated from 26 X-ray luminous galaxy clusters in Allen et al. (2004, quoted in Table 2 for 
a LCDM universe), (2) / gas (-R500) from 35 highly luminous (Lx ^ 10 45 erg s —1 ) objects in Ettori & Fabian (1999). The observational data are normalized 
to ft b h 2 = 0.0214 ± 0.0020 (Kirkman et al. 2003), H = 70 km s" 1 Mpc" 1 and Q m = 0.3. 



power-spectrum, which produces more evolved clusters (Borgani 
et al. 2005, in preparation). 

The radial distribution of the gas and baryon mass fractions 
in these massive systems present quite a similar behavior. For one 
representative object, we plot in Fig. [2] these ratios as a function 
of radius, out to 3 x -R200, at 2 — and 1, for the six physical 
schemes adopted in our analysis. The total baryon fraction mea- 
sured in the radiative models is higher than the value obtained for 
the non-radiative runs (code= G, GV) within the virial radius and 
reaches the cosmic value at about 3 x 7?2oo- The gas fraction is 
larger when the winds are stronger as a consequence of their effect 
in preventing gas removal from the hot phase. At 2 — 0, it increases 
typically with radius reaching 50 and 80 per cent of the value mea- 
sured at -R200 at r < O.I-R200 and at ~ 0.3i?2oo, respectively. At 
2 = 1, the gas fraction is less concentrated: the radii at which the 
50 and 80 per cent of / gas (< R200) is reached are a factor 1.4 and 
2 larger than at z — 0. 

The quantities Y gas and Yb at -R2500, -R500 and -R200, as a 
function of redshift and physics included in the simulations, are 
plotted in Fig. [5] In the inner cluster regions, the dissipative action 
of radiative cooling enhances the average Yb to super-cosmic val- 
ues at high redshift. At late times cooling is less efficient, and Yb 
declines, although to values (~ 0.9 at z = 0) that remain higher 
than those of the non-radiative runs. A smaller scatter and more 
widespread agreement among the different physical regimes are in- 
stead found in the outskirts (r > -R500). The gas fraction within 
i?2500 is about 0.3 times the cosmic value at z = 1 and 0.6 at 
z — 0, whereas is more tightly distributed around 0.6 — 0.7 at larger 
radii, with evidence of larger values in the presence of strong winds. 
We also compare in Fig.[3]our simulation results with the observed 
/gas distribution in highly X-ray luminous clusters at -R2500 (from 
Allen et al. 2004) and at R 500 (from Ettori & Fabian 1999). Simu- 
lations clearly indicate a sizeable underestimate of the hot baryons 
budget, both at -R2500 and at R$oo- When extra physics is added 
to the action of gravitational heating, lower hot gas fractions re- 
sult. The discrepancy with the inferred observed fraction signals 
the existence of systematic errors, either in our physical treatment, 
or in estimates of the observed fraction, or possibly both. It is worth 
noticing that total mass estimates, for instance, suffer from system- 



atic differences when measured from X-ray analysis and from dark 
matter particles in simulations, mainly owing to bias in the X-ray 
spectral temperature measurements (see, e.g., Mazzotta et al. 2004, 
Vikhlinin 2005 and, specifically related to the systematics in X- 
ray mass estimates, Rasia et al. 2005). The observed gas fraction 
measurements are well represented by the values measured in the 
non-radiative simulations also in the inner regions, where varia- 
tions by < 15 per cent due to the action of non-thermal pressure 
support (induced by the reduced viscosity scheme) encompass the 
observed variance. 

The values of, and the relation between, Y gas and Yst ar are 
shown in Fig. [4] and [5] At -R2500, the estimate of Yb is strongly 
affected by the physics involved in the accumulation of the clus- 
ter baryons, with super-cosmic values mainly due to a very high 
star-formation efficiency (Y s t a r > 0.5), particularly at high redshift 
(see Fig.|4j. The effect of this cooling excess is to sink more gas in 
the core to maintain the pressure support, thus making the overall 
baryon fraction super-cosmic and the gas fraction less concentrated 
at high redshift. At R500 and -R200, the dispersion in the estimate 
of Yb is 2-3 per cent of the measured mean (Fig.[3J. We explain the 
reduction in the scatter at large radii by the common history that 
the baryons experience in our simulations when the properties are 
averaged over a volume large enough not to be dominated from the 
physics of the core. With respect to the mean values measured at 
^2500 > Y s tar decreases by a factor > 2 and Y gas increases by more 
than 30 per cent, with peaks of 2-3 when no winds, or weak winds 
combined with the effect of conduction, affect the ICM physics 
(Fig-ID. The ratio / ga s//star, that is < 1 at -R2500 and 2 > 0.3, 
increases with both radius and redshift and becomes ~ 3 at -R200 
and z = in the presence of winds (see Fig.|5J. Once again, the 
absence of galactic winds increases the efficiency of star forma- 
tion and leaves the / ga s//star ratio below 2 within the virial radius. 
These values are well above the observed constraints, also plotted 
in Fig. [5] Lin, Mohr & Stanford (2003) compare the gas mass mea- 
sured in nearby X-ray bright galaxy clusters with the stellar masses 
evaluated from K-band luminosities of the member galaxies. From 
their estimates, converted to the total mass measurements in Ar- 
naud et al. (2005), we infer an average M sas /M sta , T w 8.7 (r.m.s. 
= 2.7) at -R500 in systems with gas temperatures larger than 3 keV, 
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that is still a factor between 2.5 and 6 larger than what obtained 
in our simulated objects (Fig. [5J, thus witnessing the presence of 
significant overcooling. 



4 COMPARISON TO PREVIOUS ANALYSES 

While the comparison between results from different non-radiative 
simulations is relatively straightforward, when extra-physics is in- 
cluded, such as radiative cooling, star formation and stellar feed- 
back, the results on the hot and cool baryons distributions is sen- 
sitive to the implemented scheme. Therefore, this difference has to 
be taken into account in performing such a comparison. 

In their non-radiative simulations, Eke, Navarro & Frenk 
(1998) found that the gas fraction within i? v ir is, on average, 87 
per cent of the cosmic value, that is reached at about 3-R v ir- For 
dynamically relaxed systems, they measure Y gaa = 0.824 ± 0.033 
at r = 0.25-Rvir « -R2500 (as quoted in Allen et al. 2004). No 
evident evolution is present between 2 = 1 and 0. A comparable 
value was obtained as average of the baryon fraction measured in 
a set of non-radiative simulations of one single cluster presented 
in the Santa Barbara Comparison Project (Frenk et al. 1999), with 
/b/(^b/^m) of 0.92 (rms: 0.07) at the virial radius. Similar re- 
sults have been obtained also in other SPH simulations (e.g. Bialek 
et al. 2001, Muanwong et al. 2002). 

These results agree with our estimates for our simulations 
of massive galaxy clusters with gravitational heating only: Y>(< 
-R200) = 0.89 ± 0.02 (0.90 ± 0.03 at i? vir ). A slight increase at 
high redshift is however measured (lb = 0.94±0.03 at z = 1). On 
the contrary, the baryon fraction decreases moving inwards, with 
Yb(< R 500 ) = 0.87 and Y b (< #2500) = 0.82 

Rravtsov, Nagai & Vikhlinin (2005) studied the baryon 
fraction in nine galaxy clusters spanning a decade in mass 
and simulated with the Eulerian adaptive mesh refinement N- 
body+gasdynamics ART code, for both non-radiative and radia- 
tive cases. For non-radiative simulations, they measure Yb ~ 1 
(at r > 3i? vir ), 0.97 ± 0.03 (at R vil ), 0.94 ± 0.03 (at i? 5 oo) 
and 0.85 ± 0.08 (at -R2500), that are consistently higher than our 
values by few per cent. This systematic difference by ~ 5% at 
r > -R2500 has been already shown by the comparison between 
the ART and GADGET codes discussed in Kravtsov et al. (2005) 
and, even though smaller than the 10% differences measured in the 
Santa Barbara cluster comparison project (Frenk et al. 1999) be- 
tween the gas fraction obtained from Eulerian and SPH codes, is 
still relevant when using simulations to calibrate systematics in the 
estimate of the baryon fraction in clusters. 

When the radiative cooling and other physical processes are 
turned on, the number of reliable comparisons that can be done is 
reduced. In the radiative simulations of Muanwong et al. (2002), 
Y gaa is between 0.6 — 0.7 in high— M systems and Yb « 0.9 at 
the virial radius. When pre-heating with an extra energy of 1.5 



keV per particle at z = 4, these authors find Y gl 
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0.8 — 0.9. Kravtsov et al. (2005) measure the gas fraction to rise 
from 0.46±0.06 (at ifesoo) to 0.65±0.06 (at R V1I ), whereas Y b re- 
mains always around 1 (from 1.22 ± 0.11 at ifesoo to 1.02 ± 0.02 
at R v ir). These values are roughly in agreement with our results, 
in particular for the set of simulated objects with weak winds: at 
2 = 0, Ygas = 0.51 ± 0.02 at #2500 and 0.70 ± 0.03 at ifeoo, 
while Yb is flat around 0.93 ± 0.02 and never super-cosmic, i.e. 
Yd < 1. Different numerical implementations of the cooling, star- 
formation and feedback processes are expected to contribute to this 
systematic difference between the predictions of ART and GADGET 



simulations. These differences are, in fact, more significant in the 
radiative runs and depend on the amount of baryons cooled in stars 
(see also Kravtsov et al. 2005). Such differences, which are nu- 
merical in origin, should definitely be considered as theoretical un- 
certainties, when using simulations to calibrate systematic biases 
in the observational estimate of the baryon fraction within clusters. 
The action of stronger winds is to increase the gas fraction, whereas 
the absence of winds reduces Y gas to 0.29 and 0.59 at -R2500 and 
R200, respectively. 



5 IMPLICATIONS FOR THE CONSTRAINTS ON 
COSMOLOGICAL PARAMETERS 

Our results have a direct implication on the systematics that affect 
the constraints on the cosmological parameters obtained through 
the cluster baryon mass fraction (e.g. White et al. 1993, Evrard 
1997, Allen et al. 2002, Ettori et al. 2003, Ettori 2003, Allen et 
al. 2004). We remind that, once a representative gas fraction, de- 
noted here / gas , is directly measured from X-ray observations and 
a statistical relation between the average / star and / gas is adopted, 
the cosmic mass density parameter can be then evaluated as 

Ybfib 



/gas (l H~ /star/ /gas) 



(l) 



where the "hat" indicates the observed quantities and the cosmic 
baryon density flb is assumed from primordial nucleosynthesis cal- 
culations or the measured anisotropies in the cosmic microwave 
background. In recent years this method has been also extended to 
the measure of the dark energy density parameter, (Oa, w, see, e.g., 
Allen et al. 2002, Ettori et al. 2003, Allen et al. 2004) under the as- 
sumption that the gas fraction remains constant in redshift (Sasaki 
1996, Pen 1997). Since the gas fraction scales with the angular di- 
ameter distance as / gas oc d A ' 5 , the best choice of cosmological 
parameters is defined as the set of values that minimizes the \ 2 
distribution of the measured gas fraction at different redshifts with 
respect to the reference value: 



r 



Yb fib^n 



f 1 + / s t ar // g 



d A (z;On = 0.3; n A = 0.7) 



(1a(z; fi m ; Qa) 



.(2) 



Despite its conceptual simplicity and straightforward applica- 
tion, this method makes some assumptions that have to be tested 
before the error bars estimated for the matter and dark-energy den- 
sity parameters can be accepted as robust and reliable determina- 
tion of both statistical and systematic uncertainties. In the present 
discussion, we highlight two of the assumptions generally adopted, 
but never verified: (1) the mean value of Yb does not evolve with 
redshift, (2) a fixed ratio between / star and / gas holds in a clus- 
ter at any radius and redshift. As we have shown here, both these 
assumptions are not valid in our simulated dataset whatever is the 
physics included in the simulations, in particular when considering 
the inner part of the clusters. Allen et al. (2004) use the simula- 
tion results by Eke et al. (1998) to fix Y b = 0.824 ± 0.033 at 
T ~ -R2500 for their sample of Chandra exposures of the largest 
relaxed clusters with redshift between 0.07 and 0.9. We notice, for 
instance, that, while this value is in agreement with our simula- 
tion results at z — in the runs with Gravitational heating only 
(Yb = 0.82 ± 0.06), it is definitely lower than what we estimate 
at higher redshift (e.g. Y b = 0.86,0.89,0.91 at z = 0.3,0.7, 1, 
respectively). This increase of Yb with redshift is the consequence 
of the different accretion pattern of shock-heated baryons at differ- 
ent epochs. At later times, accreting gas had more time to be pre- 
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Figure 4. The comparison between the outputs of the 6 different physical schemes investigated. The gas (solid line) and stellar (dashed line) mass fractions, 
normalized to the cosmic value, are plotted at i?250o(~ 0.3-f?20o)> RbOo{~ 0.7i?20o) an d J?200- Their evolution with redshift is indicated by the thickness 
of the lines (from thickest line/larger symbols to thinnest line/smaller symbols: 2 = 0, 0.7, 1). The shaded regions indicate the ltr range of y s tar, i^gas and 
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FwW FwWC 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 

r/R200 r/R 2 oo 

Figure 5. Ratios between the cumulative gas and stellar mass fractions within a given radius as a function redshift (solid line: z = 0; dashed line: z = 0.7; 
dotted line: z = 1) for the different physics included in the simulations. The large dots show the mean value (and relative standard deviation) obtained from 
the cosmological sample. The diamonds indicate the ratios expected, for a given M$oo, from the equation 10 of Lin, Mohr & Stanford (2003), based on 
near-infrared observations of massive X-ray galaxy clusters. 



shocked into filaments. As a consequence, they have a relatively 
higher entropy, thus relatively increasing the radius (in unit of the 
virial radius) where accretion shocks take place. This is the reason 
why, as shown in Fig. [2] at z = 1 /b reaches the cosmic value at 
relatevely smaller radii than at z = 0. 

Since the tighter cosmological constraints provided by the 
cluster gas fraction alone are on Q m (of the order of 16 per cent 
at la level; e.g. Allen et al. 2004), we try to quantify the effect of 
the variation of the baryonic components with the radius and the 
redshift on this estimate. To this purpose, we use equation Q and 
evaluate first how f2 m changes by varying Y b . In runs with grav- 
itational heating only, the increasing baryon fraction with redshift 
induces larger estimate of f2 m with respect to what obtained from 
local measurements of Y b '. 



S^m - H m _ AQm _ Y b (< Ra,z = Zq) _ 

On fi m Y b (<R A ,Z = 0) 

is +0.09 at Ra = -R2500 and z — 0.7 for the case with gravita- 
tional heating only and +0.11 at z = 1. (Here the prime symbol ' 
indicates the corrected value with respect to the reference one). Us- 
ing instead the runs with reduced viscosity, the deviation decreases 
to about +0.05. As for the radiative runs, the bias is of the order of 
20 per cent, that reduces to 2 per cent in presence of strong winds at 
z — 1. When outer cluster regions are mapped (i.e. r ~ -R500), the 
deviation converges to similar amounts due to the limited impact 
of cooling and feedback over large volumes: variations between 
+0.03 (weak winds) and +0.10 (strong winds) become compara- 
ble to Af2 m /f2 m ~ +0.08 as measured in non-radiative runs. 
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A further contribution to the uncertainties comes from the de- 
pendence upon the radius and redshift of the ratio /star //gas- In 
the observational determination of the baryon fraction from eqn.Q 
this quantity is generally assumed to be 0.16, as measured in the 
Coma cluster within the virial radius [e.g. White et al. 1993; note 
that our estimate from Lin et al. (2003) and adopting the total mass 
measurements in Arnaud et al. (2005) is 0.11 ± 0.04]. If we com- 
pare the ratio (f> — / s tar// g as measured for a Coma-like simulated 
cluster at -R200 and z — with the estimates at other redshifts 
(z = 0.7 and 1; see, e.g., Fig. [5}, we evaluate from equation^ 

Hm - »m = An m = l + 0(<fi A ,Z = Q) _ 1 _ _ Q Q5 _ 
fim fim 1 + </>(< Ra, Z = Zo) 

When 4> = /star //gas measured locally at -R2500 is compared 
with the corresponding value at different redshifts, the deviation 
ranges between —0.74 (when winds are excluded) and —0.38 
(when strong winds are present), whereas it is about —0.10 at 7?5oo- 

As for the runs with gravitational heating only the effect of the 
variation of Yb with redshift and overdensity implies Af2 m /f2 m < 
+0.11, thus comparable to the current statistical uncertainties from 
Chandra observations of the massive clusters out to z = 1.3 (Ettori 
et al. 2003, Allen et al. 2004). However, when the extra-physics of 
the radiative runs is included, Afi m /fi m has two contributions of 
~ +0.10 and < —0.05, due to an increase with redshift of (1) Yh 
(see Fig.|3j and (2) the stellar to gas mass fraction ratio (see Fig.[5}. 
Both these effects are caused by a more efficient star formation in 
high redshift clusters. 

In general, our results indicate that it may be dangerous to use 
simulations to calibrate observational biases for precision determi- 
nation of cosmological parameters from the gas fraction in clusters. 
Although none of our simulation models includes a fair description 
of the actual ICM physics, it is interesting that different models 
provide different redshift-dependent corrections for the estimate 
of the cosmic baryon fraction from observations of the gas and star 
density distribution within clusters. If applied to observational data, 
such corrections would induce sizeable differences in the determi- 
nation of the matter and dark energy density parameters. 



6 SUMMARY AND DISCUSSION 

We have analyzed the distribution of gas and stellar mass fraction 
in simulated massive X-ray galaxy clusters as function of (i) the 
radius expressed in the form of -R2500, Rsoo and R200, (ii) the red- 
shift at which the structure is identified, (iii) the physical processes 
determining the evolution of the baryons in the cluster potential 
wells (i.e. gravitational heating, radiative cooling, star formation, 
conduction and galactic winds powered by supernova explosions). 
Our main results can be summarized as follows. 

(i) . As for the cluster set extracted from a cosmological box, 
which are simulated by including cooling, star formation and feed- 
back with weak (340 km s" 1 ) winds, we find at -R200 Yb = 0.93, 
Vgas = 0.74 and Y st ar = 0.20, with scatter around these val- 
ues of 2, 4 and 8 per cent, respectively. These results are vir- 
tually independent of the cluster mass over the range M v ir ~ 
(0.5 — 13) x 10 14 h~ 1 Mq. The dispersion relative to the mean 
value measured at -R2500 is a factor of about 3 larger than -R200 • 

(ii) . In the four massive (M200 > 10 15 h^ 1 Mq) galaxy clus- 
ters simulated with 6 different physical schemes, we find that the 
cosmic value of the baryon fraction flb/fi m is reached at about 
3 x 7?2O0- The gas fraction increases radially, reaching 50 (80) per 
cent of the value measured at -R200 at r ~ 0.1(0.3)i?2oo at z = 0. 



At z = 1 the same values are reached at radii which are about 40 
per cent larger. This indicates that / gas tends to be less concentrated 
at higher redshift, where the more efficient star formation causes a 
more efficient removal of gas from the hot phase in the central clus- 
ter regions. We also find that in these clusters the amount of hot 
baryons, in unit of the cosmic value, is less scattered and less de- 
pendent on the particular physics adopted when it is measured over 
larger cluster regions. In the runs with Gravitational heating only, 
Yb = /bar/ (fib/fim) ranges between 0.82 ± 0.06 (at -R2500) and 
0.89 + 0.02 (ati?20o)atz = 0. It increases at z= 1 to 0.91 + 0.10, 

0. 95 + 0.04 and 0.94 + 0.03 at .R2500, -R500 and ifeoo , respectively. 
This increase of Kb with z owes to the smaller radii (in units of 
the virial radius) at higher redshift where shock-heated baryons ac- 
crete, permitting the /b values to reach the cosmic value at higher 
overdensities, for example, at z = 1 than at z — (see, e.g., 
Fig.|2j. These values are roughly in agreement with observational 
constraints obtained from highly X-ray luminous clusters (e.g. Et- 
tori & Fabian 1999, Allen et al. 2004) at -R2500 and 7? 50 o. Using an 
SPH scheme with reduced viscosity, the non-thermal pressure sup- 
port contributed by turbulent gas motions, reduces the gas fraction 
by 15 per cent at -R2500 and by ~5 per cent at r > JJsoo- 

Adding extra-physics to the action of gravity reduces the amount 
of diffuse baryons by 20-40 (20-65) per cent at R500 (-^2500). This 
is due to the high star formation efficiency that is not sufficiently 
mitigated by our assumed supernova feedback. The presence of 
strong winds (assumed here v < 480 km s _1 ) induces the fol- 
lowing effects: (1) a reduction of Ystar by about a factor of 2 with 
respect to the case with no winds; (2) a milder radial dependence 
of Ystar and Ygas ; (3) larger values for the ratio between / gas and 
/star, that is however still smaller, by about a factor 2.5, than the 
observed values (Lin et al. 2003) . 

(iii) . At i?2500, the estimate of Yb is strongly affected by the 
physics included in the simulations. Super-cosmic values are found 
for the radiative runs, as a consequence of the exceedingly effi- 
cient cooling process, which causes a rapid sinking of baryons in 
the cluster cores. Moreover, differences in the formation history of 
individual clusters contribute to a wider spread in the Yh values, 
as shown in the runs including only gravitational heating. In these 
cases, the standard deviation is of about 10 per cent. 

(iv) . At -R500 and -R200, the dispersion in the estimate of Yh is 
2-3% of the measured mean. With respect to the mean values mea- 
sured at i?2500, Ystar decreases by a factor > 2 and Yg as increases 
by more than 1.3, with peaks of 2-3 in the simulations without 
winds, or with weak winds combined with conduction at 1/3 the 
Spitzer value. 

We explain the reduction in the scatter at r > -R500 by the com- 
mon history that the baryons experience when their properties are 
averaged over a large enough volume. In the outer regions, the 
gas evolution is mainly driven by the action of gravity. The ratio 
/gas //star, that is < 1 at -R2500 and z > 0.3, increases with radius 
and redshift and becomes ~ 3 at -R200 and z — in the presence 
of winds. 

(v) . Our results have direct implications on the calibration of the 
systematic effects that limit the use of the X-ray gas mass fraction 
as a cosmological tool (e.g. Evrard 1997, Ettori et al. 2003, Allen 
et al. 2004). Owing to the uncertainties in the modelization of the 
ICM, we emphasize that our simulations cannot be used to cali- 
brate the observed / gas . However, our results emphasize the poten- 
tial problems related to the estimate of the cosmic baryon fraction 
from the gas mass fraction as measured in central cluster regions; 

1. e. the variation of Yb with redshift and a proper estimate of the 
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contribution of stars in galaxies to the total cluster baryon budget 
as a function of radius and redshift. 

For the non-radiative runs the effect of the variation of Yb with 
redshift and overdensity implies a Afi m /fi m < +0.11 that is 
comparable to the typical statistical uncertainties from Chandra 
observations of the massive clusters out to redshift z — 1.3 (Et- 
tori et al. 2003, Allen et al. 2004). However, when star formation 
and feedback are also included, Afi m /fl m has two contributions of 
~ +0.10 and < —0.05, due to an increase with redshift of (1) Yb 
(see Fig.[3} and (2) ration between stellar and gas mass (see Fig.[3}. 

Our results show that star formation in our radiative runs is 
still too efficient, especially in the cluster central regions, even in 
the presence of rather strong galactic winds. As a result, the impact 
of cooling and star formation in the simulated clusters is still much 
larger than in realistic cases. On the one hand, this calls for the 
need to introduce more effective feedback mechanisms, e.g. related 
to Active Galactic Nuclei, which are able to prevent overcooling in 
the central cluster regions (see, e.g., the effect of the injection of su- 
personic AGN jets in the ICM of simulated clusters in Zanni et al. 
2005). On the other hand, our analysis demonstrates that using sim- 
ulations to calibrate the systematic uncertainties in the estimate of 
cosmological parameters from the cluster gas mass fraction may be 
quite problematic, especially if it has to be used as a high-precision 
tool to measure the cosmic density parameters associated to matter 
and dark energy. 
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Table 1. Distribution of the baryons in the two sets of simulated galaxy clusters. The mean values of the indicated ratios are shown. 



STATUS N 2 



#2500 
^tar 



GV 



FwW 



4 0.0 0.817(0.064) 

0.3 0.864(0.057) 

0.7 0.888(0.020) 

1.0 0.910(0.094) 



FwWC 4 



0.0 
0.3 
0.7 
1.0 

0.0 
0.3 
0.7 
1.0 

0.0 
0.3 
0.7 
1.0 



0.690(0.045) 
0.720(0.069) 
0.720(0.033) 
0.722(0.106) 

0.513(0.021) 
0.448(0.040) 
0.341(0.059) 
0.288(0.036) 

0.538(0.058) 
0.448(0.070) 
0.302(0.090) 
0.247(0.052) 



0.817(0.064 
0.864(0.057 
0.888(0.020 
0.910(0.094 

0.690(0.045 
0.720(0.069 
0.720(0.033 
0.722(0.106 



0.419(0.029) 0.931(0.021 

0.565(0.054) 1.013(0.063 

0.749(0.114) 1.090(0.060 

0.884(0.098) 1.172(0.090 

0.398(0.049) 0.935(0.021 

0.552(0.072) 1.000(0.063 

0.803(0.169) 1.105(0.079 

0.965(0.160) 1.212(0.126 



#500 
^star 



0.874(0.023) 
0.907(0.022) 
0.945(0.028) 
0.947(0.037) 

0.818(0.018) 
0.865(0.019) 
0.914(0.038) 
0.912(0.045) 

0.658(0.029) 
0.664(0.025) 
0.649(0.024) 
0.631(0.036) 

0.655(0.031) 
0.652(0.026) 
0.631(0.028) 
0.620(0.058) 



0.262(0.011) 
0.286(0.022) 
0.305(0.021) 
0.330(0.038) 

0.257(0.015) 
0.281(0.017) 
0.306(0.018) 
0.327(0.048) 



0.874(0.023) 
0.907(0.022) 
0.945(0.028) 
0.947(0.037) 

0.818(0.018) 
0.865(0.019) 
0.914(0.038) 
0.912(0.045) 

0.920(0.026) 
0.950(0.027) 
0.953(0.014) 
0.961(0.028) 

0.912(0.027) 
0.932(0.025) 
0.937(0.017) 
0.946(0.025) 



#200 
^star 



0.894(0.020) 
0.926(0.019) 
0.933(0.027) 
0.942(0.033) 

0.864(0.010) 
0.900(0.013) 
0.913(0.025) 
0.935(0.016) 

0.701(0.027) 
0.708(0.018) 
0.687(0.028) 
0.684(0.029) 

0.714(0.027) 
0.726(0.026) 
0.707(0.034) 
0.704(0.033) 



0.894(0.020) 
0.926(0.019) 
0.933(0.027) 
0.942(0.033) 

0.864(0.010) 
0.900(0.013) 
0.913(0.025) 
0.935(0.016) 



0.226(0.009) 0.927(0.018) 

0.243(0.018) 0.951(0.019) 

0.259(0.013) 0.946(0.019) 

0.269(0.019) 0.953(0.033) 

0.221(0.011) 0.936(0.017) 

0.236(0.017) 0.962(0.014) 

0.250(0.008) 0.958(0.027) 

0.263(0.021) 0.967(0.027) 



F 1 0.0 


0.294(-) 


0.682(-) 


0.976(-) 


0.535(-) 


0.370(-) 


0.905(-) 


0.593(-) 


0.333(-) 


0.927(-) 


0.3 


0.169(-) 


1.083(-) 


1.252(-) 


0.491(-) 


0.447(-) 


0.938(-) 


0.594(-) 


0.369(-) 


0.963(-) 


0.7 


0.123(-) 


1.465(-) 


1.588(-) 


0.456(-) 


0.486(-) 


0.942(-) 


0.564(-) 


0.371(-) 


0.935(-) 


1.0 


0.180(-) 


1.054(-) 


1.234(-) 


0.558(-) 


0.409(-) 


0.967(-) 


0.578(-) 


0.362(-) 


0.940(-) 


FsW 1 0.0 


0.638(-) 


0.302(-) 


0.939(-) 


0.699(-) 


0.199(-) 


0.898(-) 


0.747(-) 


0.176(-) 


0.923(-) 


0.3 


0.600(-) 


0.439(-) 


1.038(-) 


0.731(-) 


0.245(-) 


0.976(-) 


0.774(-) 


0.203(-) 


0.978(-) 


0.7 


0.439(-) 


0.617(-) 


1.055(-) 


0.718(-) 


0.263(-) 


0.981(-) 


0.740(-) 


0.210(-) 


0.951(-) 


1.0 


0.409(-) 


0.552(-) 


0.960(-) 


0.745(-) 


0.248(-) 


0.993(-) 


0.730(-) 


0.207(-) 


0.937(-) 



